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Modern  radars  may  incorporate  pulse-to-pulse  carrier  frequency  modulation  to  increase 
probability  of  detection,  to  reduce  vulnerability  to  jamming,  and  to  reduce  probability  of  inter- 
ception. However,  if  coherent  processing  is  used  for  clutter  rejection,  the  frequency  of  N 
consecutive  pufses  must  be  held  constant  for  A^-pulse  clutter  cancellation  or  doppler  filtering.  If 
M pulses  are  transmitted  during  the  time  the  antenna  illuminates  a target,  there  are  M/N  coherently 
integrated  echoes  available  for  noncoherent  integration  in  the  computer  or  the  operator's  display 
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20.  Abstract  (Continued) 

to  further  improve  the  signal-to-noise  ratio  (SNR).  In  this  report,  analytical  and  simulation  methods  are 
employed  to  determine  the  balance  between  coherent  and  noncoherent  integration  that  yields  the 
greatest  SNR  improvement.  Attention  is  focused  upon  three  models  which  include  FFT  doppler 
filtering  and  different  systems  of  combining  and  noncoherently  integrating  doppler  filtered  signals. 
Curves  of  detectable  SNR  as  a function  of  M and  N are  presented  for  all  three  models. 
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FREQUENCY-AGILE  RADAR  SIGNAL  PROCESSING 


INTRODUCTION 

Modern  radars  sometimes  incorporate  pulse-to-pulse  changes  in  carrier  frequency  (a)  for 
the  sake  of  increasing  detection  probability,  (b)  to  reduce  vulnerability  to  jamming,  and  (c)  to 
reduce  probability  of  interception.  However,  coherent  processing  such  as  a coherent  moving- 
target  indicator  (MTI)  is  often  necessary  to  suppress  clutter  echoes,  which  means  that,  during 
the  period  of  time  coherent  processing  takes  place,  the  carrier  frequency  should  be  held  con- 
stant. It  is  possible,  in  principle,  to  perform  coherent  integration  of  a train  of  frequency- 
hopped  echoes.  However,  the  response  of  such  a system  to  small  changes  in  range  and 
doppler  frequency  causes  problems  beyond  the  scope  of  this  report.  Thus,  we  consider  a 
frequency-agile  radar  that  transmits  groups  of  N identical  pulses,  where  carrier  frequency  is 
fixed  within  a group  (block)  but  is  varied  from  block  to  block. 

In  this  report  we  address  the  problem  of  choosing  the  optimum  value  of  N (the  number 
of  pulses  in  a constant-frequency  block),  assuming  random  frequency  modulation  from  block 
to  block.  The  total  number  of  pulses  A/ transmitted  during  the  antenna  beam  dwell  time  is  as- 
sumed to  be  fixed.  The  method  explored  here  shows  that  a substantial  improvement  in  re- 
ceiver sensitivity  results  from  an  optimum  choice  of  /V,  especially  for  large  values  of  M. 

Fundamental  detection  theory  is  reviewed  in  the  next  section  and  is  subsequently  applied 
to  three  receiver  models,  which  are  defined  and  analyzed.  The  appendixes  provide  needed  de- 
tail to  justify  the  method  of  calculating  the  performance  of  the  receiver  models. 


PROBLEM  DEFINITION 

Here  we  consider  the  problem  of  selecting  optimum  processing  for  a sequence' of  A/ radar 
echoes,  which  are  divided  into  M/N  blocks  of  /V  pulses  per  block.  The  carrier  frequency  ,/„, 
(.see  Fig.  I)  is  a.ssumed  constant  within  a block;  random  changes  in  carrier  frequency  are  as- 
sumed from  block  to  block.  Af  and  /V  may  take  any  integer  value.  An  optimum  processor  is 
defined  as  one  that  minimizes  the  input  signal-to-noise  power  ration  required  for  given  proba- 
bilities of  detection  Pj  and  false  alarm  P^^.  If  no  further  constraints  are  added,  the  optimum 
processor  will  be  the  well-known  matched  filter  HI.  Since  the  waveform  .v(f)  is  known  exactly 
(except  for  time  of  arrival),  a linear  filter  having  an  impulse  response  /t(t)  = it.s(  —i  + Tq  ), 
would  maximize  the  signal-to-noise  ratio  at  its  output  and,  hence,  minimize  detectable  power,  S 
(the  constant,  ^ pulse  width,  and  5 is  sometimes  termed  receiver  sensitivity). 

At  present,  difficulties  in  implementation  prevent  coherent  integration  of  frequency-agile 
echoes.  Such  a processor  is  known  to  have  been  implemented  in  software  [2,3],  but  not.  to  the 
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BLOCK  1 BLOCK  2 BLOCK  M/N 


*'  *2  *N  *N»1  *2N  *M 

TIME  - 

Fig.  1 - Block-io-block  frequency-agile  radar  signal  Cf  -carrier  frequency) 

best  of  our  knowledge,  in  hardware.  Hence,  we  consider  coherent  integration  of  each  block  of 
N pulses  and  noncoherent  integration  of  M/N  blocks.  We  wish  to  find  a specific  value  of  N 
which  minimizes  the  input  signal-to-noise  ratio  ^ required  for  detection. 


PROCESSOR  MODELS 

Figures  2a  and  2b  represent  examples  of  the  processors  examined  in  this  report.  The 
processor  of  Fig.  2a  is  optimum  under  the  above  constraints,  whereas  the  second  processor 
(Fig.  2b)  has  the  advantage  of  simplicity  since  only  one  output  channel  is  required. 


SWITCHING 


THRESHOLD 


Sm  + N(TI 


REFERENCE 

SIGNAL 


(a) 


CHANNEL 

INTEGRATORS 


DETECTORS 


SITI  + NIT) 
.» 


REFERENCE 

SIGNAL 


(bl 


Fig  2 — Block'io-block  radar  signal  processing  sysiems 


A filter  matched  to  a single  doppler-shifted  block  of  A pulses  is  identified  by  an  impulse 
response  matching  the  rectangular  block  of  A pulses  (Fig.  1).  Since  the  doppler  frequency  is 
unknown,  ii  is  necessary  lo  provide  multiple  fillers  to  span  the  interval  between  doppler  ambi- 
guities. This  leads  to  the  incorporation  of  the  well-known  fast  Fourier  transform  (FFT)  as  il- 
lustrated in  Fig.  2;  both  processors  are  identical  up  through  A magnitude  detectors  operating 
on  A outputs  from  the  A-point  FFT.  It  is  assumed  that  the  filter  block,  labeled  MF,  is  matched 
to  the  waveform  of  each  pulse. 
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Deleclion  theory  ll,4|  defines  a maximuni-likelihood  detector  for  the  case  of  a signal  hav- 
ing one  unknown  parameter  to  be  a bank  of  matched  filters,  which  span  the  range  of  uncer- 
tainty (in  this  case  doppler  frequency),  followed  by  threshold  detectors  operating  on  the  filter 
output  magnitudes.  From  this  point  of  view  (maximum-likelihood  detection),  an  A'-point  FFT 
processor  is  optimum  for  detecting  a target  on  the  basis  of  a single  block  of  pulses.  We  there- 
fore incorporate,  in  the  models  of  interest  (Fig.  2),  optimum  A-pulse  coherent  processing. 

Constant-Velocity  Target  Detector 

It  can  be  argued  that  the  processor  of  Fig.  2a  is  optimum  for  the  signal  specified  in  Fig.  1. 
An  echo  from  a target  of  velocity  v will  be  integrated  in  doppler  filters  identified  by  the 
relation  = 2,/',v/c,  where  ,/j/,  is  the  doppler  frequency  corresponding  to  the  carrier  frequency 
./,  of  the  ith  block.  Since  the  echo  will  appear  in  different  doppler  channels  as  ,/,  changes, 
a switching  matrix  is  required  to  deliver  all  MIN  coherently  integrated  blocks  to  the  appro- 
priate "video"  integrator  corresponding  to  target  velocity  v.  It  is  noted  that  A F/Av  velocity 
channels  are  required,  where  A F is  the  range  of  target  velocities  and  Av  = c/2  AT,  /,  = 
X/2A7/  is  the  mean  velocity  resolution  of  the  FFT  doppler  channels. 

Although  we  do  not  prove  optimality,  this  processor  integrates  all  echo  energy  over  an 
A/-pulse  dwell  time  of  the  antenna  beam  on  a target  prior  to  combining  with  noise-only  chan- 
nels. Thus,  this  processor  is  expected  to  be  more  efficient*  than  that  of  Fig.  2b,  which  com- 
bines doppler  channels  prior  to  video  integration  in  the  peak  selector.  The  integrated  ouputs 
form  the  basis  of  an  estimate  of  target  velocity;  i.e.,  all  A//A  outputs  {|>’,  | ) due  to  a target  of 
velocity  v will  be  connected  via  the  switching  matrix  to  the  corresponding  velocity  channel  in- 
tegrator, say  v^.  Hence,  the  threshold  detectors,  set  to  select  the  largest  response,  will  most 
probably  select  the  one  containing  the  target  signal.  In  this  sense,  the  processor  of  Fig.  2a  may 
be  referred  to  as  a velocity  estimator,  where  v^  is  an  estimate  of  target  velocity  v. 

Reference  Processor 

Here  we  define  a reference  processor  which  is  identical  to  that  of  Fig.  2a  except  that  only 
one  velocity  channel  (hence,  one  video  integrator)  is  required  since  we  assume  that  target 
velocity  is  known.  The  effect  of  this  assumption  is  to  avoid  the  need  to  correct  calculated  SNR 
values  for  doppler  uncertainty.  If  the  doppler  frequency  is  unknown,  the  threshold  values 
must  be  increased,  as  A F/Av  increases,  to  maintain  a constant  false  alarm  rate. 

Peak  Filter  Output  Integrator 

The  advantage  of  processing  block-to-block  agile  signals,  using  the  system  of  Fig.  2b 
rather  than  that  of  Fig.  2a,  is  obvious;  no  switching  matrix  and  just  one  integrator  are  required. 
On  the  other  hand,  a disadvantage  stems  from  the  fact  that  no  doppler  information  is  available 
from  the  peak  output  integrator.  When  samples  of  noise  only  are  added  to  samples  containing 
signal  plus  noise,  the  signal  becomes  more  difficult  to  discern  from  noise,  and  the  increase  in 
signal  power  required  to  maintain  a given  detection  reliability  (/*j  and  Z',,,  ) is  referred  to  as 
collapsing  loss.  However,  as  will  be  shown  later,  the  performance  for  Fig.  2b  is  only  slightly 
worse  than  that  for  the  reference  processor  defined  above. 


•This  term  is  used  lo  express  efficacy  of  signal-io-noise  raiio  iroprovcmcni;  ctTiciency  in  niaihematical  siaiiMics  lakes 
on  u difffereni  meaning  noi  iniendcd  here 
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DETECTABLE  SNR  FOR  REFERENCE  PROCESSOR 

We  now  develop,  for  the  above-defined  reference  processor,  an  expression  for  the  in- 
put SNR  required  for  specified  detection  and  false  alarm  probabilities.  To  obtain  this  expres- 
sion, we  assume  that  SNR  is  improved  linearly  with  N in  the  coherent  integrator  (FFT),  and 
we  make  use  of  the  standard  curves  of  needed  (for  different  combinations  of  P,/  and  Py„) 
at  the  input  to  the  magnitude  detector  as  functions  of  / = Ml N,  the  number  of  noncoherent 
integrations.  The  calculation  of  91^.  (Pj  = 0.9,  Pj^  = 0.7  x 10“'*^)  is  plotted  in  Fig.  3,  based 
upon  the  five  statistical  models  defined  by  Swerling  (5)  for  targets  with  fluctuating  cross  sec- 
tion. Probably,  the  most  frequently  used  is  Case  I,  which  corresponds  to  Rayleigh  scan-to-scan 
fluctuations  (all  / pulses  correlated).  Cases  I and  III  differ  only  in  the  statistical  distributions 
assumed.  The  results  in  Fig.  3 clearly  illustrate  the  greater  penalty  inflicted  by  the  Case  I 
scan-to-scan  fluctuations  as  compared  with  the  Case  II  pulse-to-pulse  fluctuations  for  large 
values  of  /,  the  number  of  integrated  pulses.  Note  that  for  / = 750,  Case  II  noncoherent  in- 
tegration (NCI)  is  more  elTicient  than  Case  I NCI  by  about  8 dB. 


Fig  3 — C'ompari.son  of  coherent  and  noncoherent  integration  of 
radar  returns  from  Swerling-model  targets  (for  P,  = 0 9.  = 0.7 

X I0  '“)  (after  Rivers.  SEMCOR  Rept.  SD-7521 3-1 ) 


The  curves  of  Fig.  4 were  developed  from  Swerling  Case  II  curves  (frequency  agility 
block-to-block),  making  use  of  the  following  procedure  and  definitions.  The  detectable  SNR  at 
the  input,  is  improved  by  coherent  integration  so  that  91^  = S9t,  neglecting  losses;  hence, 
a value  of  91  may  be  found  from  Eq.  (1): 
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Fig.  4 — Required  signal-to-noise  ratio  (for  Pj  -=0.9,  =0.7  x 10"'®) 

vs  N for  fixed  values  of  M (pulses/dwelltime) 


« = (1//V)  (MIN). 

For  example,  note  that  for  Case  11  (Fig.  3),  (100/10)/10  =8—10 


(1) 

-2  JB  = di. 


DETECTABLE  SNR  FOR  PEAK  FILTER  OUTPUT  INTEGRATOR 

The  curves  of  Figs.  4 and  5,  corresponding  lo  the  processor  of  Fig.  2b,  required  a more 
complex  calculation  because  of  the  necessity  of  accounting  for  the  collapsing  loss  encountered 
when  noise-only  channels  are  combined  with  a channel  containing  both  signal  and  noise. 
Thus,  we  make  use  of  the  fact  that  a portion  of  this  processor  has  been  analyzed  for  M-ary 
communications  (4).  In  this  case,  the  block  diagram  up  through  the  peak  selector  is  identical 
lo  the  optimum  detector  for  a communications  signal  consisting  of  any  one  of  N orthogonal 
signal  waveforms.  The  probability  of  error  (i.e.,  the  probability  that  the  output  of  at  least  one 
noise-only  channel  will  exceed  the  ouput  of  the  channel  containing  the  signal)  is  computed  us- 
ing the  relation  of  Eq.  (2),  which  is  further  illustrated  by  the  graph  of  Fig.  6; 


f>,(.)  = I -/ 


+ 00  1 

r~ — 

-oo  ~J2tt 


exp 


U - (2£/A'o)'^^l^ 


JTtt 


exp  ( — v^/2 ) (/> 


,v-i 


(2) 


where  £and  AJ,  are  signal  energy  and  noise  power  spectral  density,  respectively. 


Returning  to  Figs.  3 and  4,  note  thai  the  minima  exhibited  by  the  curves  of  Fig.  4 occur 
because  of  the  change  in  slope  of  the  Case  II  curve  (solid  line)  in  Fig.  3,  as  / = ,)// A' changes 
from  small  to  large  values.  Its  slope  exceeds  that  of  the  corresponding  coherent  integration 
curve  (dashed  line)  for  small  / becoming  smaller  for  large  /.  As  A' approaches  M (Fig.  4),  / 
decreases  lo  small  values  causing  the  value  of  (Fig.  3)  to  rise  sharply,  overcoming  the  de- 
crease due  lo  greater  coherent  processing 
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Hence,  if  we  consider  a sequence  of  / output  samples  from  the  peak  selector,  on  the  aver- 
age, samples  will  be  due  to  noise  only;  clearly  (1  -/',)/ will  be  due  to  signal  plus  noise. 
Collapsing  loss  is  well  known  [6)  to  be  a function  of  the  collapsing  ratio  p,  where,  if  m is  the 
number  of  noise-only  samples  and  n is  tf  number  ol  signal-plus-noise  samples,  then 

_ m + n _ ' ~ I _ 1 /Tj 

" n ~ - F,)l  ^ ' 

Although  Marcum  17)  calculated  collapsing  loss  as  a function  of  this  ratio  p,,  his  curves  (re- 
printed by  Berkowitz  [8]),  were  computed  for  the  case  of  a nonfluctuating  target.  The  approxi- 
mation (/,,.=  ! as  discussed  in  Appendix  A,  is  therefore  used  below  to  derive  Eq.  (4). 

Since  is  a function  of  P,,  which  is  a function  of  3^  and  .V,  we  are  led  to  the  following 
equation  which  must  be  solved  by  successive  approximations: 

/?,.  (M/A) 

NS)  . «,  ^MIN)IL,  - - «• 

Equation  (4)  is  obtained  from  the  definition  of,  and  approximation  for,  collapsing  loss  (see  Ap- 
pendix A)  and  the  relationship  for  coherent  integration.  Note  that  is  defined  to 

be  a number  less  than  one. 


SIMULATION  OF  PEAK  SELECTION  PROCESSOR 

A simulation  was  implemented  to  compute  the  detection  performance  of  a processor  us- 
ing peak  selection  of  doppler  channels.  For  comparison,  the  simulation  was  also  run  for  a pro- 
cessor using  only  a single  doppler  channel.  Figure  7 is  a diagram  of  the  simulation  showing 
the  processors  and  the  signal  and  noise  generators. 


MEASURE  OF 
RMS  NOISE 
INPUT 


t I'lg  7 — Simuhiion  nl  lafsc'l  dI  ' detector  and  Miiple  dorpl>.-  v.ltanncl  detector 
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I The  processor  of  imeresi  consists  of  a doppler  filter  bank,  inii)leniented  with  an  eight- 

point  discrete  Fourier  transform  (l)FT).  Eight  input  samples  compose  a coherent  processing 
interval  (CPI)  and  produce  eight  doppler  outputs.  The  outputs  are  detected,  the  one  with  the 
largest  amplitude  is  selected,  and  the  outputs  from  lO  CPIs  are  summed  in  a noncoherent  in- 
tegrator, Finally,  the  output  of  the  integrator  is  compared  to  a fixed  threshold  and  the  detec- 
tion decision  is  made.  The  detection  decision  is  based  on  a total  dwell  time  of  M =80  samples. 
The  second  processor  is  then  implemented  as  a reference.  This  processor,  identical  except  that 
a single  doppler  output  is  used,  represents  a case  where  the  target  doppler  is  known  and  no  col- 
lapsing loss  is  suffered. 

The  inputs  to  both  processors  are  identical  and  consists  of  in-phase  and  quadrature  (/and 
Q)  components  of  signal  plus  noise.  The  signal  components  are 

1^  = A cos  (5a) 

and 

0,  = A sin  . (5b) 

where  ,/  is  the  doppler  frequency,  k the  doppler  filter  number,  and  /V' the  number  of  points  in 
the  DFT.  The  additive  noise  has  a uniform  phase  distribution  and  an  approximately  Gaussian 
amplitude  distribution  given  by 

/y  = B cos  Itt  R,  (iSa) 

C),v  = ^sin  IttR.  (6b) 

and 

12 

B = X -0-5).  (6c) 

1 

where  //  is  a uniformly  distributed  random  number  between  0.0  and  1.0. 

I The  first  step  in  the  simulation  is  to  determine  the  threshold  setting  T for  various  false 

I alarm  rates.  The  results  shown  in  Fig.  8 were  obtained  by  setting  the  input  signal  to  7ero  and 

measuring  the  probability  of  false  alarm  as  a function  of  the  threshold  setting.  Now,  we  can 
use  threshold  settings  at  a desired  and  measure  Pj  by  running  the  simulation  with  various 
‘ values  of  signal  amplitude. 

The  results  of  the  simulation  of  both  processors  are  shown  in  Fig.  *),  where  the  input 
' SNR  is  plotted  as  a function  oi  Pj  for  various  For  a =0.9  and  = 10~\  the  proces- 

sor selecting  the  largest  filter  output  requires  2.3  dB  higher  SNR  than  the  reference  processor 
For  greater  Pj  and  smaller  B,,,,  the  performance  of  the  "largest  of"  processor  approaches  that  of 
the  single  filter  case. 


C ONCLl  SIGNS  AM)  RECOMMF.NDA  I IONS 

! Wc  have  shown  that  the  choice  of  block  size  in  a block-to-block  frequency-agile  radar  i-- 

^ important  if  maximum  receiver  sensitivity  is  to  be  realized.  The  curves  of  Fig  4 exhibit  dis 

^ tinct  minima  in  signal-to-noi,se  ratio  required  for  reliable  detection,  the  locations  of  which 

1 depend  upon  the  total  number  of  pulses  available  for  coherent  and  noncoherent  processing. 
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These  data  correspond  to  a reference  case  in  which  ii  is  presumed  ihai  the  doppler  frequency 
of  ihe  larger  is  known,  so  that  collapsing  loss  is  zero.  The  collapsing  loss  corresponding  lo  the 
velociiy-esiimator  processing  (Fig.  2a'  is  easily  determined  from  the  reference  case  (to  a good 
approximation)  by  computing  the  number  of  doppler  resolution  cells  required  to  span  the 
range  of  target  velocities.  If  this  number  is  greater  than  A,  then  the  collapsing  loss  is  equal  to 
the  SNR  increase  needed  lo  reduce  by  the  factor  A'.  Here,  collapsing  loss  arises  from  the 
fact  that  more  than  one  channel  may  contribute  a false  alarm,  requiring  an  increase  in  thres- 
hold and  SNR  to  maintain  a desired  false  alarm  rale. 

Figure  4 also  contains  data  corresponding  to  peak  seleci/iniegrale  processing  (Fig.  2b)  in 
which  multiple  doppler  channels  are  combined  by  selecting  only  the  largest  output  for  subse- 
quent video  integration.  Table  1 compares  the  three  processors  for  two  values  of  total  dwell 
lime,  M = 100  and  M = 200.  where  in  each  case  the  optimum  N was  selected. 


Table  1 — Comparison  of  Processor  Performance 
(P,,  = 0.9,  P,,,  = 10  “*) 


Processor 

M 

N 

opt 

SNR  (dB) 

T , (dB) 

Reference 

100 

17 

-3.6 

0 

200 

33 

-6.8 

0 

Velocity  Estimator 

100 

17 

-2.75 

0.85 

(Fig.  2a) 

200 

33 

-5.15 

1.65 

Peak  Select/Integrate 

100 

20 

-2.1 

1.5 

(Fig.  2b) 

i 

200 

50 

-5.1 

1.7 

It  should  be  noted  that  the  SNR  penalty  for  using  the  simpler  system  involving  just  one 
integrator  is  less  than  1 dB  for  M = 100,  and  less  than  0.1  dB  for  M = 200.  However,  we 
must  consider  the  effect  of  the  approximation  used  in  calculating  performance  of  the  latter  sys- 
tem. Basically,  the  approximation  involves  two  assumptions  discussed  in  greater  detail  in  Ap- 
pendix A: 

1.  We  assumed  3?.  —3?^.  in  deriving  Eq,  (4)  used  to  compute  SNR  for  the  peak 
seleci/integrate  system. 

2.  The  effect  of  the  peak  selector  on  signal  and  noise  distributions  was  ignored;  i.e..  P^/ 

. and  Pi^  were  assumed  to  depend  only  on  the  SNR  at  the  threshold. 

I 

Table  2 compares  the  data  calculated  from  Eqs.  (I)  and  (4)  with  Monte  Carlo  simulations 
of  the  reference  and  the  peak  select/inlegrale  processors.  No  approximations  are  involved  in 
Eq.  (1)  except  those  used  by  Marcum,  and  the  signal  in  the  development  of  Table  2 is 
nonffuctuaiing.  The  results  clearly  indicate  that  the  simulation  and  the  calculalion  for  the  refc; 
cnce  processor  agree  within  0.2  dB;  the  simulation  for  the  peak  selector  is  higher  by  0.5  dB 
than  the  c^tputed  result  involving  = 1 — P^.  and  1.1  dB  for  the  case  invoKing 
Z.J  = xA*—  Pr-  These  collapsing  loss  expressions,  described  in  Appendix  .A.  hold  for  largs 
and  small  SNR.  respectively.  It  might  be  inferred  from  Table  2 that  the  effect  of  the  second 
assumption  (mentioned  above)  is  lo  yield  optimistic  results;  hence,  C,  = 1 — P^  provides  a 
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more  accurate  estimate  of  SNR.  This  would  tend  to  agree  with  the  expected  influence  of  the 
peak  selector  upon  the  noise  statistics;  i.e.,  it  is  reasonable  to  expect  the  distribution  change  to 
increase  A*,,;  It  has  been  suggested  that  L^.  depends  more  heavily  on  A/ than  on  M,  the  total 
number  of  integrated  pulses  (9). 

Table  2 — Comparison  of  Performance  Estimates  Obtained  from 
Eqs.  (1)  and  (4)  and  Corresponding  Monte  Carlo  Simulation* 


Processor 

Method 

1 

SNR  (dB) 

Reference 

Eq.  (1) 

-4.9 

Reference 

Simulation 

-^.1 

Peak  Select/Integrate 

Eq.  (4)  (L  =1  - P ) 

^ c r 

-3.3 

Peak  Select/Integrate 

Eq.  (4)  (Z.^,  = Vl  -P,) 

-3.9 

Peak  Select/Integrate 

Simulation 

-2.8 

•P^-0.9,  P - to  *,  M - 80,  iV  - 8 

a fa 


In  conclusion,  it  should  also  be  noted  that  the  processor  of  Fig.  2a  provides  an  estimate 
of  target  doppler  as  well  as  nearly  optimum  sensitivity.  Further  attention  should  be  given  to 
this  technique.  For  example,  a measure  of  reliability  of  velocity  data  for  targets  of  various 
speeds  and  accelerations  could  be  defined  and  tested.  Although  the  utility  of  frequency-agile 
systems  depends  on  the  range  of  frequencies  available  and  interference  with  other  users  of  mi- 
crowave bands,  doppler  ambiguities  are  reduced  and  blind  speeds  are  eliminated  by  this  tech- 
nique, so  that  an  estimate  of  radial  velocity  is  available  on  each  scan  of  the  antenna.  Velocity 
data  should  tend  to  reduce  track-while-scan  errors  and  the  hand-off  delay  to  weapons  systems. 
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Appendix  A 

COLLAPSING  LOSS  ASSOCIATED  WITH  PEAK  OUTPUT 
SELECTION  AND  INTEGRATION 

In  this  appendix,  we  are  interested  in  finding  the  value  of  ^ needed  to  provide  the  re- 
quired detection  reliability  measured  by  and  (detection  and  false  alarm  probabilities)  for 
the  system  illustrated  in  Fig.  Al.  Since  N-1  noise-only  channels  are  combined  with  the  chan- 
nel containing  the  signal  in  a peak  selector,  the  SNR  delivered  to  the  integrator  is  degraded;  to 
compensate,  must  be  increased  by  an  amount  defined  here  to  be  the  collapsing  loss  L^.. 
The  necessary  increase  in  may  be  approximated  by  introducing  the  diagram  of  Fig.  A2. 
Here,^^  is  the  SNR  needed  If  all  integrated  pulses  contain  signal  plus  noise  and  we  define  col- 
lapsing loss  as 

L,  (Al) 

Note  that  the  output  SNR,  is  the  same  for  both  diagrams  (Figs.  Al  and  A2). 
Although  it  is  known  that,  in  general,  equal  detection  reliabilities  (F^,  Pfa)  will  not  be 
achieved  by  making .^g  the  same  in  Figs.  Al  and  A2  because  of  the  differing  statistical  distri- 
butions caused  by  the  effect  of  the  peak  selector  (Fig.  Al),  comparisons  of  data  computed  from 
the  results  of  the  following  analysis  indicate  that  this  effect  is  small. 


DEMODULATORS 


Fig.  Al  — Noncoherent  integration  of  peak-selected  doppler  filtered  signals 
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We  assume^ihat  SNR  improves  linearly  with  /,  because  of  video  integration,  so  that,  in 
f'ig.  A signal  is  assumed  to  be  present  in  just  one  of  the  .V  doppler  chan- 

I -I 

nels  and  SNR  is  assumed  to  be  Equal  amounts  of  noise  are  assumed  in  all  channels.  On 
the  average,  integrated  pulses  are  due  to  noise  only  (i.e.,J?.,  = 0)  where  is  the  so-called 
error  probability  described  in  the  foregoing  sections  and  in  Appendix  C.  Hence,  if  ^^y  is  as- 
sumed constant  for  / pulses,  then  = ,5?.  for  the  remaining  (1  — terms,  and 

(I  -P,)i 

^0  = Z (A2) 

I -I  ‘ 

In  Fig.  A2,^^.  =S^qII  = (1  — ).?^.,  since ^0  same  (Eq.  (A2))  as  in  Fig.  A1  and 

is  assumed  to  be  constant. 


Now,  to  evaluate  =^yl^y,  we  introduce  an  expression  connecting  with  Si^. 
Equation  (A3),  given  by  Davenport  and  Root,*  makes  clear  the  nonlinear  dependence  of  on 
The  same  expression  connects^,  and  ^y. 


= kp 


1 +20ly  ' 


(A3) 


where  =1  for  a rectangular  sine  wave  signal  envelope.  Equation  (A3)  holds  for  Gaussian 
noise,  a square  law  detector,  and  where !%.  is  defined  as  the  ratio  of  output  signal  and  noise- 
only  powers.  Although  the  equation  for  5?,  varies  with  assumptions  regarding  detector  law  and 
SNR  definition,  we  are  dealing  here  with  a mechanism  affected  only  slightly  by  the  functional 
dependence  of  on  ^y.  In  the  following  paragraphs,  this  is  demonstrated  by  comparing 
values  of  obtained  by  calculations  based  upon  two  assumptions:  (a)  — ^y  and  (b)  as 

given  by  Eq.  (A3). 


Let  = where  an  example  of  the  function  / is  given  in  Eq.  (A3).  Note  that 

the  collapsing  loss  may  be  expressed  in  terms  of  the  inverse  function  /"'.  Since 
^y  = / “*  ),  we  may  substitute  this  expression  into  Eq.  (Al),  and  we  have 

/■'(^j  /-'  ((1  -P,)^-] 

L, ; ^ ; (A4) 

/ (.^;  ) / (•'^.  ) 

where  the  right-hand  side  of  Eq.  (A4)  makes  use  of  the  substitution  = (1 
derived  from  Eq.  (A2). 


If  we  let  Pfi.  ~ fi^.y  then  Eq.  (A4)  reduces  to  the  form  introduced  in  the  text  in  Eq.  (4); 
i.e.,  Lj.  =1  —P^.  If  we  require  greater  accuracy  (especially  for  .'3?,  < < 1/2),  we  must  use 
Eq.  (A5)  obtained  from  Eqs.  (A2)  and  (A4). 


i9ly  {\  -Py) 

1 + 2. My 


1 -f 


(A5) 


*W,  B.  Davcnporl,  Jr.  and  W L Root,  An  Inirmliniion  in  the  Thmn  nl  Rumlnm  SikiiciR  nnj  Smsc.  Mctiraw -Hill,  New 
York,  1958,  p.  266. 
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In  solving  Eq.  (A3)  two  roots  result;  the  root  leading  to  negative  values  of  was  dis- 

carded. 

Equation  (AS)  is  plotted  in  Fig.  A3  for  several  values  of  P^.  Note  that  the  approximation 
Vl  ~f‘r  — holds  for  small  values  oi  < < Ml  with  L,.  = I - for  ^y  > > 1/2. 


Fig.  A3  — Collapsing  loss  (Eq.  (AS))  vs.  signal-to-noise  ratio 
and  error  probability  P 


Appendix  B 

SIGNAL-TO-NOISE  RATIO  CALCULATIONS 
FOR  SWERLING  CASE  II 


Original  works  by  Marcum*  and  Swerlingf  have  been  discussed  by  many  authors.  A re- 
port by  L.  F.  Fehlner§  includes  calculations  for  the  pulse-to-pulse  fluctuation  case,  Swerling 
Case  II  (SW-2).  All  five  of  their  statistical  models  are  illustrated  in  the  curves  of  Fig.  Bl.  The 
purpose  of  this  appendix  is  to  present  expressions  for  Pj  and  for  SW-2. 


NUMBER  OF  PULSES  INTEGRATED.  I 

Fig.  Bl  — Comparison  of  coherent  and  noncoherent  in- 
tegration of  radar  returns  from  Swerling-model  targets  (for 
P.  -0.9,  P,  —0.7  X lO'”*)  (after  Rivers.  SFMCOR  Kept 
SD-75213-1) 


*J  1 Marcum,  "A  Statistical  Theory  of  Target  Detection  by  Pulsed  Radar,  Rand  Research  Memo  RM-754,  Dec  1947, 
and  RM-653,  July  1948  (reprinted  in  IRE  Trans.  IT-6,  (No.  2),  59-267  (1960). 

tP  Swerling,  Probability  of  Detection  for  Fluctuating  Targets,"  Rand  Corp  Research  Memo  RM-I2I7.  Mar  1954  (re- 
printed in  IRE  Trans.  IT-6,  269-308  (Apr.  1960). 

§L.  F.  Fehiner,  "Target  Detection  by  a Pulsed  Radar,"  Johns  Hopkins  Applied  Physics  Laboratory  Report  TO  451,  Jul 
1962 
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Appendix  C 

COLLAPSING  RATIO  CALCULATION 


Van  Trees*  discusses  "M-ary"  detection  of  signals  in  additive,  white  Gaussian  noise.  The 
expression  is  presented  in  Eq.  (Cl)  (also  in  Eq.  (2)  in  the  text)  for  error  probability  /*,((). 
This  quantity,  used  in  this  report  to  compute  the  collapsing  ratio  and  collapsing  loss  L^,  is 
given  as 


/>,(*)  -1  - 


dx 


exp 


U-  (2£/yVo)'^^l^ 


/: 


exp  ( -y^iDdy 


N-\ 


(Cl) 


The  double  integral  in  Eq.  (Cl)  has  been  evaluated  extensively  by  Urbano.t  Tables  of 
values  from  his  report  were  utilized  in  our  analysis  as  follows: 


Urbano’s  tabulated  function  (a)  is  identical  to  the  second  term  of  (Cl)  if  substitutions 
for  Vand  a are  made  and  if  the  variable  of  integration  is  changed  appropriately.  Thus, 


M -1 

dy\  dz. 


Letting  z ^ x — a results  in 

By  inspection  of  (Cl)  and  (C3),  we  may  write. 


\f- 


-y^n 


yfYn 


S-\ 


dx. 


(C2) 


(C3) 


P,U)  -1  -/’/v(a)  (C4) 

if  we  let  a = -JIE/Nq,  where  £ is  the  energy  in  each  input  signal  waveform  and  Vq  is  the  in- 
put noise  power  density. 


*H  L.  Van  Trees.  Deiecnon,  Estimation,  and  Modulation  Theory,  John  Wiley  and  Sons,  Inc.,  New  York,  1968 
t R tt.  Urbano,  "Analysis  and  Tabulation  of  the  M-Positions  Experiment  Integral  and  Related  Error  Function  In- 
tegrals," AFCRC  Report  TR-55-100,  Apr.  1955 
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